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We report a new multinuclei based nuclear magnetic reso- 
nance method which aUows us to image the staggered polar- 
ization induced by nonmagnetic Li impurities in underdoped 
(Oe.e) and slightly overdoped (O7) YBa2Cu306+H above Tc- 
The spatial extension of the polarization ^imp approximately 
follows a Curie law, increasing up to six lattice constants at 
80/S' at Os.e in the pseudogap regime. Near optimal doping, 
the staggered magnetization has the same shape, with ^imp 
reduced by a factor 2. ^imp is argued to reveal the intrinsic 
magnetic correlation length of the pure system. It is found to 
display a smooth evolution through the pseudogap regime. 



I- INTRODUCTION 

In high Tc cuprates, strong antiferromagnetic (AF) 
correlations reminiscent of the AF ordered phase were of- 
ten argued to explain the anomalous normal state proper- 
ties of the Cu02 planes and to lead to the d-wave pairing 
in the superconducting state [1]. Nuclear magnetic res- 
onance (NMR) and neutron scattering experiments both 
showed the existence of such correlations and their de- 
crease with increasing doping. However, no agreement 
has been reached about quantitative determinations of 
the AF correlation length ^ values and about a possi- 
ble incommensurability. In underdoped YBa2Cu306-i-j;, 
some neutron experiments detect a magnetic response 
peak with a temperature (T) independent q width, which 
corresponds to ^ « 2.5 (in lattice units) [2]. But this 
response could be dominated by a macroscopic inhomo- 
geneity in the large single crystals used or according to 
other experiments an incommensurability, so that ^ could 
be underestimated. In contrast, the analysis of Pines et 
al (noted hereafter MMP) of the Cu NMR relaxation 
times Ti and T2 suggests a Curie- Weiss behavior for 
^ which might reach values as large as 7 at lOOiC [3]. But, 
various reexaminations of the MMP analysis lead to sub- 
stantial uncertainties in their ^ values, so that at Og.Oj 
a model-free determination of ^ versus T is still lack- 
ing. Near optimal doping, the magnetic response peak 
measured by neutrons broadens a lot and its intensity is 
reduced indicating smaller magnetic correlations. This 
makes the experimental data difficult to analyse quanti- 
tatively. 

We propose to use nonmagnetic impurities to deter- 
mine this length scale ^ in cuprates. Nonmagnetic sub- 
stitutions on the magnetic sites of low dimensionnal AF 
correlated systems are known to lead to the appearance of 



an AF staggered magnetization. This has been observed 
by NMR in cuprates [4] [5] [6], spin chains [7], and spin 
ladders [8]. In the case of an AF S = \ spin chain, the in- 
duced polarization could be imaged site by site [7] . Its ex- 
tension was found to be identical to the correlation length 
of the pure system. In the same spirit, we propose a de- 
tailed quantitative study of the staggered magnetization 
induced by nonmagnetic Li and Zn substituted at the Cu 
sites of Cu02 planes in YBa2Cu306+y. Previous studies 
proved that the staggered polarization has a non-Fermi 
liquid T dependence, confirming the correlated nature of 
the planes [5]. However, no quantitative determination of 
the extent or shape of this polarization could be achieved. 
It is not yet clear whether the decay of the spin polariza- 
tion can be described by a single length scale. One can 
indeeed wonder whether charge and spin degrees of free- 
dom can result in different responses. Here, we measure 
the modifications induced by nonmagnetic impurities on 
the NMR spectra of all the available nuclei, ^^O, and 
the ''Li impurity itself. Combining all these informations 
we determine, for the first time, the shape of the impurity 
induced polarization on a significantly large lengthscale 
and deduce that a single lengthscale £,imp can be used. 

This "multinuclei" analysis allows us to deduce the 
spatial variation of the polarization with temperature 
and doping. We demonstrate that S,imp increases with 
decreasing temperature in the pseudogap phase. This po- 
larization persists, but is reduced in O7, out of the pseu- 
dogap phase, together with the Kondo-type T variation 
of its magnitude evidenced in [9] . Comparison with the 
MMP analysis suggests that ^imp is, as in spin chains, 
a direct measure of the correlation length ^ of the pure 
system. 

The paper is organized as follows: Sec. II is devoted to 
the description of samples and of the NMR experimen- 
tal conditions. In Sec. Ill we decribe the principle of 
our multinuclei method. In Sec. IV, we present the ex- 
perimental results. The next two sections are concerned 
with the description and analysis of the experimental re- 
sults. We detail in Sec. V the steps for determining the 
induced magnetization and present our results for the in- 
duced magnetism at Og.e and O7. A comparison with 
neutron experiments in Zn substituted O7 is discussed 
in Sec. VI. We compare also quantitatively our findings 
with the neutron scattering and nuclear spin relaxation 
measurements in the pure YBCO compound. 

II- EXPERIMENTAL DETAILS 

NMR measurements presented hereafter were done on 
YBa2(Cui_a;„Li2;„)306+j, compounds with xn = 0, 0.01 
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and 0.02. According to measurements done in [9] on the 
same samples, a quantity Xpiane = 0.85 x a; at of Li is found 
to substitute the Cu site of the planes. For all the sam- 
ples, the ^^O enrichment procedure was carried out as 
in [5]. The samples were fully reoxidized at 300°C lead- 
ing to the slightly overdoped composition y = 1 (noted 
hereafter as O7). Deoxidation of part of each batch al- 
lowed to produce the underdoped samples ?y = 0.6 (Oe.e) 
under thermogravimetric control. The powders were fi- 
nally aligned in Stycast epoxy in a high magnetic field 
so that their c axis could then be chosen to lie parallel 
to the NMR static field. Our samples present a misori- 
entation of less than 10° of the c axis which introduces 
a second order quadrupolar broadening of the Cu spec- 
tra but has negligible effect on the ^^Y spectrum and on 
the central transition of the ^''O NMR. We performed 
i7q^ 89y^ 7-|^j measurements at O7 and Oe.e (some of the 
and ""Li NMR data at Oe.e were already published 
in [9]). Spectra were obtained using Fourier Transform 
of the echoes obtained with a standard (7r/2) — t — tt 
sequence under typical fields Hext of 7 Tesla. We also 
performed ^^Cu NMR measurements in the underdoped 
compound at 300K and lOOK in a 5Tesla applied field. 

Ill- MULTINUCLEI METHOD 

The key point of our multinuclei method is that each 
nuclear species is coupled diff(;rently to the electronic spin 
density on the Cu sites, i.e., they have different hyper- 
fine form factors [10]. This was already judiciously used 
in the pure system to distinguish between uniform and 
antiferromagnetic spin fluctuations through Ti measure- 
ments [3] . We proceed in the same spirit for the impurity 
problem. We assume that an impurity positioned at the 
origin induces a polarization {S^{~r',T)) noted 5™ on 
the Cu site at = (m, n) in lattice units from the im- 
purity. The "'^^O nucleus is predominantly coupled to the 
spin density of its two copper neighbors and has an NMR 
resonance shifted by : 

''K{m,n)='' a(5r + 5-+i) 

where Ahfg/Hext, ^''^^hf is the hyperfine cou- 

pling between ^^O and one Cu site in the c direction of 
the measurements and g is the Lande factor. The other 
nuclei follow in the same manner: 

^^K{m, n) = ^'>a{S:^ + 5™+! + 5™ ^ + 5™+/), 

63x(m, n) = "^M^Si: + -t- + + 

'^K{0, 0) = * 4 * 

where A and B are respectively the Cu on-site and near- 
est neighbor supertransferred hyperfine couplings [11]. 
Due to its specific position, the ''Li shift probes only the 
polarization of the r = 1 Cu near neighbors (noted nn). 
For the other nuclei, the different geometrical factors act 
like filters. To illustrate this effect, in Fig.l we plot the 



magnitude of these shifts | K{m,0) \ along the (1,0) di- 
rection and on checkerboards for an exponentially decay- 
ing polarization. For an alternated polarization, the ^''O 

shift ^'^K{m, n) measures a quantity which roughly corre- 
sponds to the \st derivative of the polarization envelope. 
Similarly, the ^^Y shift ^^K{m, n) rather represents the 
curvature of the polarization envelope. The more com- 
plex filtering operated by the shift measurement of ''"'Cu 
combines both a direct contribution of the spin polar- 
ization and of its derivatives. From these couplings, it 
can be intuitively understood that the various nuclei will 
not be sensitive to the same range of the polarization. 
This can be seen in Fig. 1, where, for a given polariza- 
tion, we represent the local shifts using the proper filter 
for each nucleus. Furthermore, for each nuclear species, 
the existence of various soiirces of NMR line broadening 
in the pure materials limits the experimental sensitiv- 
ity and defines a cut-off on the distance which can be 
probed. This leads us to conclude that ^^Y primarily 
probes the nearfield (r < 3) of the polarization whereas 
^''O is sensitive to r < 8. The comparison between NMR 
spectra of these different nuclei should then allow to re- 
cover an information on a large range of distances. We 
note that ^"^Cu NMR in principle probes the full range 
and might thus appear the ideal probe to determine the 
whole polarization. Unfortunately the experimental con- 
ditions are not sufficient to exploit the ^"^Cu NMR spectra 
for a quantitative analysis, as will be detailed hereafter. 

IV- EXPERIMENTAL RESULTS 

Our measurements show indeed qualitatively different 

features for these various nuclei. This confirms that these 
nuclei do not probe the same spatial ranges of the Cu spin 
polarization. In underdoped Og.e samples, two well sepa- 
rated satellites appear in the *^Y NMR spectra as shown 
in Fig. 2. They correspond to the first and second nn to 
the Li as anticipated from Fig. 1. The ■'■''O NMR spectra 
do not show such satellites, but are symetrically broad- 
ened by the impurities as shown in the inset of Fig. 3. 
This broadening has a faster T— variation than that of the 
*^Y satellites shift which follows a Curie-like variation. 
The ''Li NMR shift also follows a 1/T Curie law in the 
T— range measured here [9]. As for the ^^Cu NMR , the 
large second order quadrupole effects makes it extremely 
sensitive in pure samples on the preparation conditions 
and on the quality of the alignement of the powder. So 
it is hard to compare accurately the spectra of pure and 
substituted samples. We have however measured between 
300i4r and lOOif an absolute increase of the linewidth 
of 23000 ± 4000ppm for 2% Li. Although these numbers 
will be useful for comparison with the data on the other 
nuclei, we considered that a systematic study of the Cu 
NMR in our specific experimental conditions would not 
be useful. This will be discussed further in Sec. V B 1. 

In slightly overdoped O7, the impurity effects on the 
NMR spectra are not only strongly reduced but also dis- 
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play a qualitative change in their T dependences. Dis- 
crete NMR satellite resonances are no more detectable 
in the ^^Y spectra which only exhibit a slight broadening 
which varies roughly as l/T (see the inset of Fig. 2). The 
^'^O width also increases as l/T (see Fig. 3) while the 
^Li shift follows a l/(r+135) Curie- Weiss behavior. The 
behaviors of the ^^Y and the ^^O spectra both at Oe.e 
and O7 are qualitatively very similar to those observed 
in Zn-substitutcd compounds. Their magnitude only dif- 
fer by at most 20% from that observed on Li substituted 
samples for the same in-plane impurity concentration. 
Therefore, the observations done for Li impurities can be 
considered as universal signatures of nonmagnetic impu- 
rity effects in cuprates. 

V- DETERMINATION OF THE INDUCED MAG- 
NETIZATION 

We now apply the multinuclei analysis to these findings 
in order to reconstruct the induced polarization. We per- 
formed extensive numerical simulations in order to fit the 
data, constrain the shape of the corresponding polariza- 
tion, and deduce its amplitude and the T dependence of 
its extension. Our methodology is detailed in Sec. V A, 
in the regime where the measurements show the largest 
impurity-induced changes, i.e., in the Oe.e compound at 
the lowest temperatures. We demonstrate that a single 
length scale and shape actually describes the spatial de- 
pendence of the polarization. In Sec. V B this methodol- 
ogy has been systematically applied versus temperature. 
It allows us to present our main results, i.e., the spa- 
tial characterization of the induced magnetism from an 
underdoped compound Oe.e up to a slightly overdoped 
compound O7. 

A. Methodology 

To get quantitative information from these experimen- 
tal observations, let us write the local AF polarization as 
such: 

{S^{-r,T)) = -Sl{-ir+^fT{r) (1) 

where /t("^) is the envelope of the staggered po- 
larization at a temperature T which is normalized to 
/t(T^) = 1 for II 11= 1. 

The hypcrfinc couplings of the pure compounds are 
known for all nuclei [11] and can be retained for nuclear 
sites which are not nearest neighbors of the impurities. 
However, the hyperfine coupling ^'^A\^ between the ^^Cu 

and ^^Y nearest neighbors of the same Li impurity is 
expected to be modified. The ''Li hyperfine coupling Ah / 
is not knwown accurately as well. The analysis done 
below will allow us to constrain their values. 

We have performed numerical simulations of the ^^Y 
and ^^O spectra for various choices of amplitude S'q and 
shape /t(^) on a checkerboard of 150 * 150 cells, with a 
concentration Xpiane of randomly distributed impurities. 



In Sec. V A 1, we demonstrate that the experimen- 
tal features of the ^^O and ^'^Y NMR can only be ac- 
counted for by a continuous exponentially shaped polar- 
ization and hence by only one length scale. In Sec. V A 
2 the ^''O spectrum shape further allows us to determine 
Sq, and therefore ^^Aj^f and '^Ahf. 

1. Shape and continuity of the spin polarization 

At a given temperature Sg is initially unknown and the 
^^Y first satellite NMR shift is not useful here to analyse 

the polarization shape. At Og.e? the final amplitude S'q 
and shape /t(~^) must account for (i^ the presence of 
two ^'^Y satellites on the spectnmi, the shift of the 
second ^^Y satellite and (iiij the Lorentzian like shape of 
the ^^O line. 

We shall illustrate on the %QK data the analysis done 
at various temperatures. First we use the experimental 
features of the ^®Y spectrum to constrain the polariza- 
tion envelope at short distances from the impurity. The 
assumption of a given amplitude S}^ at SQK leads to a 
unique choice for fri'r') for which only two ^®Y satel- 
lites distinct from the central line exist. This can be 
clearly seen in Fig. 4. For a given amplitude Sq, dif- 
ferent curvatures A, B and C lead to ^^Y spectra with 
qualitatively different spectral structures. Only one of 
them {B) presents two satellites and fit the second satel- 
lite shift. The decay C is too smooth to produce two 
satellites distinct from the central line. The sharper cur- 
vature A leads to a third satellite not observed experi- 
mentally. Moreover, this polarization shape cannot ac- 
count quantitatively for the second satellite position. On 
the contrary, the exponentially decaying B accounts well 
for both the presence of two satellites and the value of 
the shift of the second satellite. Hence *^Y data allows to 
constrain /t("^) up to the fourth Cu sites (r < -s/S) for a 
given choice of Sq. Anyhow, different solutions for differ- 
ent S'q are still compatible with the experiments. This is 
exemplified in Fig. 4 where both B and D polarizations 
represented in the inset of Fig. 4 lead to ^^Y spectra 
with only two satellites and the second one at the proper 
position. 

In order to constrain the shape fri'r') at further dis- 
tances, we use now the ^''O spectra. Let us consider as 
an example the short range polarization envelope rep- 
resented by the curve B in Fig. 4. For r > V5) we 
free the shape in order to account for the experimen- 
tal spectrum. Figure 5 represents simulated ^^O spectra 
and their corresponding polarizations. The only satis- 
fying fit is obtained for an exponential curvature which 
prolongates the short distance decay B. A polarization 
with two different length scales results in a more com- 
plex spectral shape. This is exemplified by the simulated 
spectra in Fig. 5. The experimental smooth lorentzian 
spectrum can only be explained by a continuous exponen- 
tially shaped envelope. On the spectra corresponding to 
discontinuous polarizations, a satellite or shoulders ap- 
pear on the ^''O spectra . This allows then to exclude 
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the scenario in which a large induced moment appears 
on a short length scale and yields a longer range spin 
polarization. 

2. Amplitude 

To determine now the polarization amplitude Sq, we 
attempt fits of the experimental ^^O spectrum with the 
various possible short range polarizations which corre- 
spond to distinct Sq such as D and E (see Fig. 6). We let 
the polarization shape free for r > v^. The Lorentzian 
^^O experimental spectrum can be fit at best by the rep- 
resented curvatures. The shape required for r > \/5 is 
again found exponential and continuous with the short 
range curvature. But all these cases do not fit the ^''O 
spectrum equally well. Case D results in an -'^^O spec- 
trum with shoulders and case E can not explain quan- 
titatively the broadening. Eventhough D, B and E all 
account for ^^Y spectra only B accounts as well for the 
^^O shape and width at the same temperature. 

Thus I'^O NMR allows to select the proper amplitude 
Sq, and therefore yields a value for the Li hyperfine 
coupling "^Ahf through Sl{T = 80) =^ K{T = 80K) 
Hext/4:'^ Afijg. Taking into account the experimental ac- 
curacy, we found that "^A^f ranges from 0.85 to 1.05kOe. 
This result is consistent with our previous upper exper- 
imental estimation "^Ahj < 2.5kOe [9], done from the 
comparison with the measured susceptibility. The value 
Sq also allows us to determine the first satellite shift 
and the first nn ^^Y hyperfine field ^^Alj = 0.8 x^^Ahf, 
which only slightly differs from that of the bulk. 

For increasing T . the ^^Y experimental spectrum al- 
ways presents two satellites and the ^^O lineshape re- 
mains Lorentzian. Hence the polarization envelope shape 
has to be exponential at all temperatures, and a single 
length scale ^i„ip is necessary to describe it. In the fol- 
lowing, we choose the Bessel functional form frir) = 
KqIt / $_imp) / Kq{1 / ^imp) , which approximates an expo- 
nential exp{—2r / S^imp) for the experimental range of val- 
ues of r/^imp explored. This mathematical notation is 
used here rather than the exponential one as it is exactly 
that which appears in the MMP analysis of the relax- 
ation data. The fact that the same analytical definition 
is taken for ^imp and ^ in the two approaches will allow 
us to compare them quantitatively. 

B. Results 



As shown above the experimental data constrains a 
Bessel-type polarization shape at each temperature. Its 
extension S^imp is then directly determined by the ^''O 
impurity-induced broadening. 

1. Induced magnetism in underdoped O^.e 



The corresponding findings for $,irnp{T) are presented 
in Fig. 7, the upper (lower) bounds of the error bars 
being associated with the higher (lower) limit for "^A^f. 
Our central results consists then in the restricted range 
of a couple of values for the amplitude Sq and the ex- 
tension ^imp{T) of the staggered polarization. From the 
obtained values of ^imp(r), we could compute the ^^Y 
satellites positions over the full range of temperatures. 
Those are found to agree reasonably well with the ex- 
perimental ones as can be seen in Fig. 2. This confirms 
a posteriori the validity of our analysis and of the value 
obtained for ^^^jjj • 

To check further the coherence of the analysis we com- 
puted the corresponding ^^Cu spectra and compared 
them to the experimental ones. At lOOK, our simula- 
tions lead to a ^^Cu broadening between 22000ppm and 
25000ppm. This value is very similar to our measured 
broadening. But a quantitative comparison between the 
experimental Cu spectra and the simulated ones requires 
to pay further attention to some experimental compli- 
cations. Recent Nuclear Quadrupole Resonance experi- 
ments evidenced that nuclei up to at least the fourth nn of 
the impurities are undetectable due to relaxation effects 
or to an impurity contribution to the electric field gradi- 
ent [12]. As such effects should then occur as well in our 
NMR measurements, we measured the xjv dependence 
of the ^^Cu NMR intensity. Assuming that a nucleus 
either fully participates to the NMR line, or is merely 
wiped out, we estimated a cutoff radius of about four 
unit cells (slightly smaller than that of [12]). We then 
computed the expected impurity-induced broadening of 
the Cu spectrum when this wipe out is taken into ac- 
count. This broadening is reduced to a value between 
16000 and ISOOOppm, slightly smaller than the experi- 
mental value of 23000ppm. Other ^^Cu measurements 
have been done on a single crystal sample with a similar 
Zn concentration but with slightly larger hole doping [6] . 
They show an impurity-induced broadening of 12600ppm 
at lOOK, almost 2 times smaller than ours. For this data 
taken in different experimental conditions (external field, 
sample doping, and dopant), comparisons beyond order 
of magnitude agreement cannot be performed. 

To conclude this comparison with ^'^Cu NMR data an 
overall agreement is obtained, but the experimental com- 
plications mentioned above prevents getting any further 
information from the NMR spectra in our experimental 
conditions. 

2. Induced magnetism in slightly overdoped Or 

As the hyperfine couplings are usually found indepen- 
dent of doping in cuprates, the value of "^Ahf determined 
from the Og.e analysis can be kept, so the data for 
K{T) directly yield the values of Sq. Hereagain, the 
Lorentzian shape of the ^^O NMR line constrains an ex- 
ponential shape for the polarization envelope, and the 
impurity-induced broadening directly allowed us to de- 
duce ^impiT), which is plotted in Fig. 7. The corre- 
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spending ^^Y simulated satellite positions are plotted in 
the inset of Fig. 2. They appear much less shifted from 
the main line than at Og.e- This is expected, as the Li 
shift indicates that the polarization on the first nn Cu 
is much smaller than in Oe.e- A convolution with the 
typical lineshape of the pure sample allows us to confirm 
that such satellites cannot be resolved in agreement with 
the data, but contribute to the broadening of the line 
which fits the experimental one. 

3. Comparison with the macroscopic susceptibility 

data 

Finally we compute the impurity contribution to the 
macroscopic susceptibility Ximp which is given by the 
sum of (S'^(T^)) on all Cu sites of the checkerboard. At 
06.6) the computed Xmip follows a Curie- Weiss law in 
piffl [T + 6) with 9 < IQQK and an effective magnetic 
moment pe// = 1.5(3)^b per impurity. This behavior 
agrees within our uncertainties with the Curie law mea- 
sured in Zn-substituted high purity compounds [13]. At 
O7, the calculated Ximp is also compatible with these 
data, as it presents a strong reduction mostly due to 
an increase of Q. This confirms the robustness of our 
analysis. This variation of the Curie -Weiss temperature 
detected on the first nn has been attributed to a Kondo- 
type screening of the induced moment by the carriers [9] . 
Prom Eq. (1), we find here that this Kondo-type T de- 
pendence merely multiplies the whole AF response. This 
situation is similar to the standard Kondo effect in dilute 
alloys where (5-(r^,T)) - (^4^p„,.^,(^))./(r) [14]. How- 
ever in our case the impurity itself is not magnetic and 
the function /(r) actually depends on T as a consequence 
of the variation of ^. 

VI- DISCUSSION 



A. Comparison to other impurity studies 

Let us compare now our findings to the existing ex- 
perimental studies. Our induced magnetization model is 
compatible with previous NMR measurements [4] [5] [6] 
in YBCO with nonmagnetic Zn impurities. Inelastic neu- 
trons scattering (INS) measurements performed on O7 + 
Zn have evidenced that Zn substitutions enhance low en- 
ergy AF fluctuations [15]. The q width of the inelastic 
peak at the AF wave vector Qaf reveals a length scale 
^„ of 1.3 cell units at IhOK and 1.6 cell units just above 
Tc- This length scale obtained either with a Lorentzian 
or Gaussian fit of the peak is very similar to our result 
(Fig. 8). The neutron peak results from the superpo- 
sition of the impurity induced AF fluctuations on the 
contribution of the pure system. It is then difficult to 
attribute a simple significance to On the contrary, 
in our analysis, we deal only with the impurity-induced 
magnetism. The quantitative agreement between these 



results is therefore quite surprising but might indicate 
that these length scales are not that different, at least 
for O7. On the other side, the longitudinal spin lattice 
relaxation Ti of ""Li [16] relates directly to the low fre- 
quency response probed by INS. The present determina- 
tion of Ahf allows us to estimate the electronic fluc- 
tuation time T of the nn induced magnetic moment. It 
corresponds to an energy scale h/r = (18.4±3.7)mey at 
lOOif in O7 (independently of Li concentration) [17] very 
similar to the value 15 ± 3meV observed by INS in O7 
with 1.6% Zn [18]. Thus at O7, in presence of 1.6% Zn, 
the fluctuation time of the induced moments is also coher- 
ent with the energy of the impurity-induced AF fluctua- 
tions measured by INS. We however stress here that the 
present NMR experiments allow to separate unambigu- 
ously the impurity-induced magnetism from the pure one. 
Our method gives the spatial dependence of the staggered 
magnetization and its temperature dependence, and al- 
lowed us to probe as well the induced magnetism in the 
underdoped regime. 

Various theoretical studies show that the AF corre- 
lations in the pure compound are responsible for the 
appearance of an AF polarization near a spin vacancy. 
Weak coupling AF spin fluctuation models lead to a po- 
larization similar to ours within RPA [19]. Strong cou- 
pling models such as t-J or Resonating Valence Bond 
also anticipate the observed polarization, at least quali- 
tatively [20]. However, contrary to the case of undoped 
spin chains none of these approaches have attempted to 
establish a relation between the extension of the induced 
magnetism and the correlation length in the pure com- 
pound. 

Let us now compare the polarization extension 

S,imp near the impurities to the correlation length ^ as ob- 
tained from experimental studies on the pure compounds. 

B. Comparison to the MMP emalysis 

Indirect information can be obtained through the anal- 
ysis of Ti and T2 ^^Cu NMR measurements as done by 
MMP. At O7, the MMP analysis yields a nearly constant 
^ ~ 2 between lOOK and 2Q0K, represented on figure 8. 
Within the error bars, this value is very close to our value 
of ^imp- In contrast, in the underdoped regime, we find a 
value of ^imp at room T about 5 times smaller than that 
deduced by MMP for ^, and a much larger T dependence 
(see figure 8). However, the MMP analysis relies on T2 
which can be strongly affected by the occurence of the 
intrabilayer coupling and also by a possible incommen- 
surability. These effects both lead to an overestimate of 
^, especially at 06.6 [21] [22]. Hence, the only reliable 
estimates from the MMP analysis are those done for O7. 
At this doping, the similarity between ^mmp and our 
S,imp is therefore a good indication that the impurity re- 
veals the correlation length of the pure compound. The 
much better experimental accuracy obtained for ^imp (T) 
as compared to £,mm p implies that ^ nearly varies as 1 /T 
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at O7. At Oe.ej our determination of S,imp{T) should then 
presumably also represent the T variation of ^. 

C. Comparison to INS measurements in the 
pure compounds 

INS experiments in pure compounds in principle al- 
low to probe ^ through the q width of the AF peak 
((1/Ag) oc when detectable. Incommensurabilities 
of the AF response might yield a broadening which leads 
then to an overestimate of Aq. Hence neutrons usually 
give a lower bound for ^. Our results are compatible with 
previous INS results if taking into account a possible in- 
commensurability or energy dependance of the neutrons 
data (see Fig. 8, in [2]). But the absence of any T varia- 
tion of the AF peak seen by INS a priori contrasts with 
our observations. 

First, these INS experiments probe high energy fluctu- 
ations, contrary to NMR which measures the zero-energy 
limit. Aeppli et al. in pure LSCO (almost optimally 
doped [23]) showed the existence of an energy dependence 
for Aq. When the energy is lowered, Aq is reduced at low 
T and enhanced at high T. This suggests that the corre- 
lation length measured at low energy could still behave 
similarly to our data for ^imp- A recent theoretical work 
[24] proposes a spin excitation function x" (g, ui) with a 
small plateau around Qaf which would reconcile as well 
the T-independent q width of the neutrons peak in the 
pure system with our T-dependent £,imp- 

Of course an incommensurability 6 of the magnetic 
response, as evidenced in pure underdoped compounds 
by some neutrons experiments [25] might as well explain 
this apparent discrepancy. We considered whether such 
a value for 6 woiild play a role in the analysis of our re- 
sults. This should yield a large-range modulation of /t(?') 
which should vanish at a distance of about eight or nine 
unit cells. Fortunately, this effect does not affect the stag- 
gered magnetism at short distances probed by ^^Y and 
^^O. We found that the changes induced in our simulated 
^^O spectra would not be detected within experimental 
accuracy. So our determination of {S^{r, T)) for r < 8 is 
therefore unaffected. In contrast ^^Cu nuclei probe the 
polarization at larger distances than ^^Y and ^^O and 
might be sensitive to such an incommensurability. But 
probing the staggered magnetization at r > 8 might only 
be ensured after securing that the observations are not 
affected by interactions between impurities. Determining 
the shape of the tail of the polarization with ^^Cu NMR 
thus requires a much better experimental resolution, i.e., 
high quality single cristals with x plane << !%• 

In summary, we have compared the impurity induced 
polarization extension (,imp{T) to the correlation length 
^(T) estimated in the pure compounds. We found that 
these quantities are similar to that obtained from the 
MMP analysis in the optimally doped samples. Various 
possible sources for a broadening of the INS AF peak in 
the underdoped samples have been proposed. This lets 



us suggest that the spatial range Cimp(T) of the induced 
magnetization around an impurity might be the best es- 
timate of ^(T). 

VII-CONCLUSION 



In conclusion, we have used NMR to determine the 
shape and magnitude of the staggered polarization in- 
duced by a nonmagnetic impurity. We have found that 
the decay of the polarization is characterized by a single 
length scale, and is best fitted by an exponential. This de- 
tailed information obtained gives a new basic constraint 
on any microscopic model of cupratcs. This experiment 
allows us to demonstrate that the AF correlations in the 
normal state of YBCO persist up to optimal doping with 
no qualitative abrupt change when crossing the pseudo- 
gap line. The extension ^j^p of the staggered polarization 
is found sizeable and T dependent in both the under- 
doped and slightly ovcrdopcd regimes. We suggest that 
this quantity might be the best estimate of the actual 
correlation length in the pure systems, although further 
theoretical support for this proposal is still required. 
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FIGURE 1 : Magnitude of the local shifts | K(m, n) \ 
at the 1^0 (b), ^^Y (c) and ^^Cu (d) nuclei for H//c from 
an impurity-induced polarization S'™ as represented in 
absolute values in (a) as function of the distance to the 
impurity along the (1,0) direction. The checkerboards 
(insets) show the same quantities on the bidimensionnal 
Cu02 layer. Black color is for the maximum intensity 
and each square corresponds to one nucleus. In (b), the 
^^O neighbor sites to the impurity (hatched area) - not 
observed up to now - are not considered here as they only 
sense one Cu site, and should undergo a large shift. 

FIGURE 2 : Data for the shift with respect to the 
main ^^Y line of the two satellites resonances detected at 
Oe.e- A reference spectrum taken at IQQK is displayed. 
The computed satellite shifts with the parameters values 
given in Fig. 6 are shown as gray lines. The inset displays 
the equivalent figure for the O7 case, and shows that the 
satellites cannot be resolved. 

FIGURE 3 : Broadening of the NMR spectrum of the 
planar oxygens induced by Li impurities with xjv = 2% 
for Oe.e (triangles) and O7 (circles). The broadening of 
the line in pure compound has been substractcd. In the 
inset is displayed a typical Lorentzian spectrum for Oe.e 
with xn = 2% Li at 80K. The small "bump" in the low 
frequency tail corresponds to the resonance line of the 
apical site, almost completely attenuated by dynamical 
contrast (fast repetition time). 
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FIGURE 4 : Comparison of the experimental ^ Y spec- 
trum at 80K in Og.e to simulated ones and their corre- 
sponding polarization (inset). Dotted Kne indicates the 
experimental satellite shift. Polarizations A, B and C 
have the same amplitude at r = 1 but different curva- 
tures. They lead to the computed spectra A, _B, and C. 
For this given amplitude, only the B exponential decay 
account for the experimental features. For a different 
r ~ \ amplitude, the polarization D leads to simulated 
spectrum D which can also account for the experimental 
shape. On these simulated spectra, the first satellite po- 
sition does not coincide with the experimental, as we use 
the value ^'^Ahf of the pure compound for the hyperfine 
coupling between and ^^Cu nn of the Li. 

FIGURE 5 : Comparison of the experimental ^'^O spec- 
trum (thick gray hue) at 80if in Og.e with Xn = 2% of 
Li to the simulated ones with the corresponding polar- 
izations in insert, with same symbols. The short range 
polarization is fixed to the B case of Fig. 4. The shape of 
the polarization for r > \/5 is varied to fit the experimen- 
tal spectrum. A convolution of the computed spectra by 
a Lorentzian function corresponding to the broadening of 
the pure sample has been performed. 

FIGURE 6 : Comparison of the experimental ^^O spec- 
trum at 80if in Oe.e with xm = 2% of Li to simulated 
ones (b) with the corresponding polarizations in (a), with 
the same symbols. For these different amplitudes D, B 
and E, the best fit of the experimental spectrum is ob- 
tained with an exponential curvature. Among these three 
cases, only B can account both for the shape and the 
broadening of the spectrum. 

FIGURE 7 : Polarization extension £^imp (as defined in 
the text) as function of temperature for underdoped Oe.e 
(closed circles) and slightly ovcrdopcd O7 (open circles). 
The inset shows the polarization envelope at SQK for Og.e 
and O7. 

FIGURE 8 : Comparison of our result for the polar- 
ization extension £,imp to the correlation length ^ in pure 
compounds extracted from INS experiments ( [2]) and the 
MMP analysis ( [3]). We plot also at O7 the characteris- 
tic length ^„ of the additional AF fluctuations observed 
in presence of Zn on the INS spectra ( [15]). 
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